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, In quantum cosmology the closed universe can spontaneously nucleate out of 

the state with no classical space and time. The semiclassical tunneling nucleation 
probability can be estimated as P ~ exp(— a 2 /A) where a=const and A is the 
^ \ cosmological constant. 

In classical cosmology with varying speed of light c(t) (VSL) it is possible to 
solve the horizon problem, the flatness problem and the A-problem if c = sa n with 
s=const and n < —2. We show that in VSL quantum cosmology with n < — 2 the 
semiclassical tunneling nucleation probability is P ~ exp(— (3 2 K k ) with /3=const and 
k > 0. Thus, the semiclassical tunneling nucleation probability in VSL quantum 
cosmology is very different from this one in quantum cosmology with c=const. In 
particular, this one is strongly suppressed for large values of A. 



1 Introduction. 

Among the most interesting alternatives to inflation in physics nowadays are truly the 
cosmological models with varying speed of light (VSL) [I], [2] 1 . In simplest case the speed 
of light c = c(t) varies as some power of the expansion scale factor: c(t) = sa n (t), where 
constant s > 0. Summarizing some of the promising positive features of these models: 

1. It can solve the horizon problem if n < 0. 

2. It can solve the flatness problem in a radiation-dominated early universe if n < —1. 

3. In case of n < — 2 the VSL models can solve the A-problem in a radiation-dominated 
early universe while inflation models can't handle it without the aid of the anthropic 
principle. 

Of course, these VSL models result in some shortcomings and unusual (unphysical?) 
features as well Pj: 

1. It is not clear how to solve the isotropy problem. 

2. The quantum wavelengths of massive particle states and the radii of primordial black 
holes can grow sufficiently fast to exceed the scale of the particle horizon. 

3. The entropy problem: Entropy can decrease with increasing time. 

1 ln fact, there are many articles about this matter. But we'd like to restrict ourselves to consider only 
those ones which has been used in this work. 



Keeping in mind all the above-mentioned problems we'd like, nevertheless, to consider 
VSL quantum cosmology. There are known the three ways to describe quantum cosmol- 
ogy: the Hartle-Hawking wave function jlj, the Linde wave function [Hj, and the tunneling 
wave function [0]. In the last case the universe can tunnel through the potential barrier 
to the regime of unbounded expansion with semiclassical probability P ~ exp(— a 2 /A). 
For the universe filled with radiation this result can be obtained in framework of a stan- 
dard WKB wave function for a particle in a potential U(a) = c 2 a 2 (l — Aa 2 /3), with the 
small energy E, see Fig.l. We have two Lorentzian regions (0 < a < a-, a > a^) and 
one Euclidean region (a[ < a < a^). The second turning point a = a,i corresponds to 
the beginning of our universe. If A = then U (a) has the form of parabola and we 
get only one Lorentzian region (see Fig.l). In this case, the universe can start at a = 0, 
expand to a maximum radius and then recollapse. If E — > then the single Lorentzian 
region is contract to the point. This, of course, comes to an agreement with the tunnelling 
nucleation probability: P — > as A — ► 0. 

In this article, however, we'll show that in quantum cosmological VSL models the 
situation can be opposite, viz: the probability to find the finite universe short after 
it's tunneling through the potential barrier is P ~ exp(— /3(n)A a ( n )) with a(n) > and 
/3{n) > when n < —2 or for — 1 < n < —2/3. After the tunneling one get the finite 
universe with "initial" value of scale factor a« ~ A -1 / 2 , so the probability to find the 
universe with large value of A and small value of a, is strongly suppressed. The reason 
of this is that, for the case A — > 0, the potential U(a) is transformed into the hyperbola 
which is located under the abscissa axis and thus such a universe can start at a ~ the 
regime of unbounded expansion (see Fig. 2). As a result, we get the single Lorentzian 
region which is not contract to the point at E — > 0. 



2 Albrecht-Magueijo-Barrow VSL model 

Lets start with the Friedmann and Raychaudhuri system of equations with k — +1 (we 
assume that G=const): 
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where a = a(t) is the expansion scale factor of the Friedmann metric, p is the fluid pres- 
sure, p is the fluid density, k is the curvature parameter, A is the cosmological constant, 
Co is the modern value of speed of light (3 x 10 10 sm/sec) and ao is the modern value of 
the scale factor (ao ~ 10 28 sm). 
Using (Jlj) one get 
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Choosing w — 1/3 one can solve (j2J to receive 
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where M > is a constant characterizing the amount of radiation. It is clear from the 
(pTJ) that the flatness problem can be solved in a radiation-dominated early universe by 
an interval of VSL evolution if n < —1, whereas the problem of A-term can be solved only 
if n < —2. The evolution equation for the scale factor a (the second equation in system 
|H)) can be written as 

p 2 + U(a) — E, E > 0, (4) 
where p = —aa is the momentum conjugate to a, E = 8nGM/3 and 
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The potential © has one maximum at a = a e = a/3/ (2A) such that 

s 2 3 n+1 

U e = U(a e ) = 2 n + l A n + l( n+1 )( n + 2 )' ^ 

so U e > if (i) n < — 2 or (ii) n > — 1. The first case allows us to solve the flatness 
and "Lambda" problems. The surplus dividend of the model is the presence of finite time 
region during which universe has accelerated expansion (see Appendix). 



3 The semiclassical tunneling probability in VSL mod- 
els with n < — 2: the case E < [/ e 

One can choose n = — 2 — m with m > 0. Such a substitution gives us the potential (JHJ 
in the form 
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The equation (jlj) is quite similar to equation for the particle of energy E that is moving 
in potential (JZJ), hence the universe in quantum cosmology can start at a ~ 0, expand 
to a maximum radius and then tunnel through the potential barrier to the regime of 
unbounded expansion with "initial" value a = a^, see Fig. 2. The semiclassical tunneling 
probability can be estimated as 



P~ exp ( --2 / | i>{a) r/(, j . (8) 

with 

|p(a)| = |^b(a)|, |p(a) |= ^(a) - E <U e 

It is convenient to write E = U e sin 2 with < < 7r/2. 



For the case E <C t/ e one can choose 



l/m 



and evaluate the integral (JHJ) as 
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The integral (fTT|) can be calculated for the m E Z . For example 

x /3 

ii(0) = ( 3 + cos2 ^) ( 191 ~ 78 cos2 9 + 15 cos4 ) ^6 + 2 cos 2 ~ 0.148 + 0(9 6 ), 

h ~ 0.025, I 3 ~ 0.007, J 4 ~ 0.002, 

and so on. One can further show that I m (9) >0at0<#<Cl. Thus, it is easy to see 
from (fT6|) that the semiclassical tunneling probability P — > for large values of A > 
and P -> 1 at A -» 0. 

Note, that the case c=const can be obtained by substitution m = — 2 into the (fTKll . 
Not surprisingly, this case will get us the well known result P ~ exp(— 1/A) (see |2j). 

In quantum cosmology the closed universe can spontaneously nucleate out of a state 
with no classical space and time. It mean that we must choose 

«) 2 < a 2 Pl (t) = (12) 

Substituting (JHJ) into (jT^J) one get the inequality, 

where Co is the modern value of speed of light and a the modern value of scale factor (or 
the size of visible universe: ~ 10 28 sm). For example 

Ai = 0.14xl0- 90 sm" 2 , A 2 = 0.48xl(T 80 sm" 2 , A 3 = 0.22xl0~ 74 sm" 2 , A 10 = 10~ 63 sm" 
For the case m<Clwe get 
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soA« 0.16 x 10~ 116 sm" 2 . 

It is convenient to introduce the cosmological parameter Q —\ p A (t ) \ /p c (t ), where 
p A {to) is the contribution of A to modern value of density and p c (^o) = 10 -29 gramme/sm 3 
is the critical density. One get 
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Q 1 = 0.2 x 10~ 34 , Q 2 = 0.5 x 10" 24 , fi 10 = 0.6 x 10~ 7 , Q 2 o = 1.3 x 10" , 



and so on. Hence, in framework of this model, the whole contribution of A in modern era 
is a quite negligible quantity (And that's just how it should be.). 



The obtained values of A are seemingly unnatural. Of course, we can assume that 
a[ > dpi, in spite of the (fT2j) . In this case, we receive the following picture: there is the pre- 
Big-Bang universe which can start its evolution at a ~ a Pi , expand to a maximum radius 
a£ and then tunnel through the potential barrier to the regime of unbounded expansion 
i.e. turn into the universe we see nowadays. The probability of this process will be large 
for small value of A. But the A-term in our universe will be the same as in pre-Big-Bang 
universe. Thus, the above scenario doesn't explain the reason why the A was so small 
in pre-Big-Bang universe and how the pre-Big-Bang universe was born. Therefore, if we 
want to describe the quantum nucleation of universe we must use the (JT2T) . As we have 
seen, the requirement for this is the verity of (fT3|) . i.e. the smallness of A and it is highly 
uncommon that for such values of A, the probability of quantum nucleation, in fact, will 
be so large. 

4 The semiclassical tunneling probability with n < 2 

and n > — 1 

In the case of general position the semiclassical tunneling probability with n = — 2 — m 
has the form 
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where 



F m (z, 9) = -2 m+1 sin 2 9 z 2{m+1) + 2 x 3 m (m + l)z 2 - m3 m+1 , (15) 

z is dimensionless quantity and the turning points, i.e. two real positive solutions 

of the equation F m (z, 6) = for the given 9 (it is easy to see that the equation F m (z, 6) = 
does have two such solutions at < 9 < ir/2). 

If m is the natural number then the expression (|T4| has a more simple form. For 
example 
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P 1 ~ exp - — — _ / —J( z *-z?)(z?-z*) , (16) 
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This expression can be calculated exactly: 
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with 



T . Q . 1 /2sin(0/2)\ 5 r8A 4 -13A 2 + 8 TT / 2 7T. . \ . 4 x2 . 2 
m = W5 [-^f- ) [ CQs2(g/2) n -\ arcsm ,) - 2 (2A 4 - A 2 + 2) K (/ , 2 



where /i 2 = cos9/ cos 2 (6*/2), A = cot(6 l /2), II and K are complete elliptic integral of the 
first and the third kinds correspondingly [8j. 



Similarly, 



exp r^y/( z2 + z N z2 > ( 18 ) 



18nGh L, z 



where 



9 tt\ ,39 \ 3 9 ix 

sin 6/ V 3 6 / V sin 9 3 \ sin 9 \ 3 6 



and so on. 

Therefore the probability to obtain (via quantum tunneling through the potential 
barrier) the universe in the regime of unbounded expansion is strongly suppressed for 
large values of A and small values of the initial scale factor Oj = \/3/(2 sin(6'/2)v / A). In 
other words, overwhelming majority of universes which are nascent via quantum tunneling 
through the potential barrier © have large initial scale factor and small value of A. 
Furthermore, it's clear that this case permits to obtain the analog of the (JT3l) too. 

Now, let us consider the case (ii), when n > — 1. The "quantum potential" has the 
form 

"<«> = sV ™ " s^ij) • < 19 > 

where m = n + 1 > 0, see Fig. 3. The points of intersection with the abscissa axis a are 
a = and a\ = ^/3(m + l)/2Am. Choosing E = in equation (J3J) and substituting (|T9|) 
into the (JHJ) we get 
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with Zi = A/3(m + l)/2m 2 . Thus, we have the same effect as if < m < 1/3. 
The density p is 

M 3(1 -m)s 2 s 2 (l-m)A , . 
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The second term is the curvature term, whereas the third term is the A-term. It is easy 
to see that at large a the first term falls off faster than the curvature term which in turn 
falls off faster than the A-term. At the beginning, when t = t\ and a(ti) = a\ we have 



p{ti) = -p\(h 
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2 The starting value z — means that the Universe tunneled from "nothing" to a closed universe of a 
finite radius a\ — Zx/yh, 



In other words, we have accelerated expansion at the beginning. As times goes by, the 
A-term will be dominating term so 

(, \ l/(l-m) 

with « = 1/(1 — to) G (1, 1.5), at t — ► oo. Therefore for £ — > oo we have a < 0. It means 
that accelerated expansion continue during a finite time. 



5 Peculiar cases with n = — 1 and n = — 2 

Now, lets consider the cases of n — — 1 and n = —2. The formula (|Tijl is not valid in 
these cases (to = — 1 and to = 0) so we shall consider these models separately. 
If n — — 1 (to = —1) then 

_ M As 2 3s 2 a 

therefore 

^M£>-^ + 1 > (22) 

where a* are constant and [a»]=sm. The potential (1221) has one maximum at a = o e = 
a/3/(2A) such that U e = U(a e ) = 2s 2 log(a e /a*), so if a e > a* then i7 e > (see Fig. 4) 
. We choose a* = A _1/2 . This gives us U e = 0.41s 2 > 0. For the case E < C/ e the 
semiclassical tunneling nucleation probability is 
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where the turning points are = 0.721, Zi = 1.812. As we can see from the (J2"HJ) . when 
n = — 1 we receive the aforementioned effect again. 
If n = — 2 (to = 0) then 
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We choose a* = l/(a\/A), where a is a dimensionless quantity. Thus 

C/(a) = -s 2 + ^ log (aaVX) + Q . (24) 
The maximum of potential (f24l) is located at the same point a e and 
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Therefore, C/ e > if a < 2e" 3 / 4 / v / 6 ~ 0.386. Choosing a = 0.286 and E < C/ e gets us 
the turning points ^ ~ 0.77 and ^ ~ 2.391. The potential is pictured on the Fig. 5. 
At last, the semiclassical tunneling nucleation probability is 
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6 Conclusion 



As we have shown, the semiclassical tunneling nucleation probability in VSL quantum 
cosmology is quit different from the one in quantum cosmology with c=const. In the first 
case this probability is strongly suppressed for large values of A whereas in the second 
case it is strongly suppressed for small values of A. Nevertheless, we can't really say that 
VSL quantum cosmology results in solution of the A-mystery. The problem is the validity 
the WKB wave function. And what is more, we have omitted all preexponential factors 
which can be essential ones near the turning points. 

But, all in all, the difference between P in VSL and usual quantum cosmology seems 
very interesting and very significant. 
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Appendix 



Lets take n 
JJJ) yields 



-2 — m, for the m > 0. Substituting |2*T| into the first equation of system 
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Thus we have the following situation: 
1. If 

2 3m(m + 2) 



< af <C 



then the curvature term is the dominating one and a < 0. 
2. If 
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then the dominating term is A-term and a > during this time. 
3. If 
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then the radiation term is the dominating one and a < 0. 

There are two way to interpret the region (A2). The first way is to conclude that we 
have cosmological inflation in early universe. This is possible when 0<m<l. In this 



case we can evaluate the number of e-foldings AiV during the region (A2) as 
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log m ~ — 2mAiV, A » — ~ ^ . (A3) 

If AiV ~ 60 then m ~ 0.029; if AiV ~ 100 then m ~ 0.0197. To evaluate £ one can use 
the well-known expression for the Friedmann integrals [H], 
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Since A(l/3) = 3 x 10 36 sm 2 /sec, we get E = 0.9 x 10 73 sm 4 /sec 2 . The substitution of 
A(l/3) into the (A3) results in 



A > 0.435 x 10~ 61 sm" 2 , A > 0.296 x 10~ 61 sm" 2 , 



for the AiV = 60 and AiV = 100. 

But do we really need inflation in the VSL models? The question is not quite clear. 
On the one hand, VSL models can solve fundamental cosmological problems (horizon and 
flatness problems) without inflation - and what is more, these models can solve A-problem 
whereas inflations can't do it without the anthropic principle. On the other hand, the 
simplest case of VSL cosmological models, which is the subject of this article, is facing 
with the isotropy problem (BJ. But, as we have seen, VSL model results in inflation with 
exit naturally so it will be incorrectly to oppose VSL models and the inflation. 

Another way to interpret the region (A2) is to identify this region with modern accel- 
eration of universe. This is possible if m is sufficiently large. Let us make a crude guess. 
According to modern observations we can write do/ao = 5.6nGp c /3 where p c = 10~ 29 
gramme/sm 3 . If the modern value of ao ~ a (see the inequality (A2)) then 

a 3mE 

A = -— —. (A4) 

From the (Al) we have 

2(m + l)As 2 

if the A-term is dominating one. Substituting (A4) gets us 



Ac 2 , 0.35m 4 / 3E 



8nGp c m + 1 V 5.6iiGp c 

Thus, if m > 1 then Q A ~ 0.35, A ~ 0.2 x 10" 56 sm" 2 , a ~ 10 27 sm. And these values 
are seems to be quite reasonable. 
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